Ion acceleration in multispecies targets driven by intense laser radiation pressure 
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The acceleration of ions from ultra-thin foils has been investigated using 250 TW, sub-ps laser 
pulses, focused on target at intensities up to 3 x 10 20 W cm -2 . The ion spectra show the appearance 
of narrow band features for proton and Carbon peaked at higher energy (in the 5-10 MeV/nucleon 
range) and with significantly higher flux than previously reported. The spectral features, and their 
scaling with laser and target parameters, provide evidence of a multispecies scenario of Radiation 
Pressure Acceleration in the Light Sail mode, as confirmed by analytical estimates and 2D Particle In 
Cell simulations. The scaling indicates that monoenergetic peaks with more than 100 MeV/nucleon 
energies are obtainable with moderate improvements of the target and laser characteristics, which 
are within reach of ongoing technical developments. 

PACS numbers: 



Significant attention has been paid lately to laser- 
driven ion acceleration, which potentially offer a com- 
pact, cost-effective alternative to conventional sources, 
with potential applications in scientific, technological and 
healthcare applications [l|. Most experimental research 
so far has dealt with the Target Normal Sheath Accel- 
eration (TNSA) mechanism [l|, where ions are acceler- 
ated by space charge fields set up by relativistic electrons 
at the target surfaces. TNSA ion beams have typically 
broadband spectrum, modest conversion efficiency and 
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large divergence, and reported scaling trends E oc I 
[l|, (where E and Iq represent maximum proton energy 
and peak laser intensity) which however may not apply to 
higher intensities regimes Q. A different mechanism, Ra- 
diation Pressure Acceleration (RPA) [3-16] is currently 
attracting a substantial amount of experimental and the- 
oretical attention due to the predicted superior scaling 
in terms of ion energy and laser-ion conversion efficiency. 
In this prospective, the so called 'Light Sail'(LS) [8l4l6| 
scheme, where, in a sufficiently thin foil, the whole laser- 
irradiated area is detached and pushed forward by the 
Radiation Pressure, is particularly promising. 

Among the attractive features of the LS mechanism 
are a favorable scaling with the laser fluence d, [l3| , 
natively narrow energy bandwidth, a reduced divergence, 
and a similarly efficient acceleration for both protons and 
higher mass ions, as predicted by numerous computa- 
tional and analytical studies. Experimental evidence of 
RPA-LS is however scarce, as TNSA dominates in stan- 
dard interaction conditions. Henig et. al [16] have re- 
ported on spectral features, which have been associated 
to RPA effects - namely a polarization-dependent modi- 
fication of the carbon ion spectral profile at moderate ion 



energies. 

This Letter presents experimental evidence of narrow 
band spectral features emerging from thin foil irradia- 
tion by sub-PW laser pulses. In particular, Carbon ion 
peaks of up to ~ 7 MeV/nucleon (cut-off energy in ex- 
cess of 10 MeV/nucleon) are produced for the first time, 
with nearly an order of magnitude higher particle flux 
than previously reported [16l4l8j. The spectral features, 
and their scaling with the laser and target parameters, 
point to a multispecies scenario of LS acceleration as de- 
scribed in [III . The possibility of achieving spectral 
peaks beyond 100 MeV/nucleon, a key requirement for 
hadron therapy [20[ , by tuning currently achievable laser 
and target parameters is discussed on the basis of the ex- 
perimental scaling, and supported by 2D Particle in Cell 
(PIC) simulations. 

The experiment was carried out employing the PW 
arm of the VULCAN laser system at the Rutherford Ap- 
pleton Laboratory, STFC, UK. A schematic of the ex- 
perimental setup is shown in the Fig. QJa). The laser de- 
livered ~200 J (varying from 175 J to 225 J in different 
shots) energy on target in pulses of 700 fs - 900 fs FWHM 
duration after being reflected off a Plasma Mirror (PM), 
resulting in intensity contrast ratio of 10 9 between the 
main pulse and the ns long amplified spontaneous emis- 
sion (ASE). The laser was focussed on targets at nor- 
mal incidence angle by an //3 off axis parabolic mirror. 
The laser peak intensity on the target was varied from 
7 = 5 x 10 19 W cm -2 to 3 x 10 20 W cm -2 by varying 
the focal spot size of the laser on target. The polarisa- 
tion of the laser on the target was varied from Linear (LP) 
(e = 0), to elliptical (EP) (e ~ 0.47 ± 0.02) and nearly 
circular (CP) (s ~ 1.14 ±0.04) by employing a zero order 
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FIG. 1: (a) Schematic of the experimental setup, (b) Ion spectra obtained from 100 nm Cu target irradiated by a LP laser 
pulse at io 3 x 10 20 W cm -2 . Different line colors correspond to different ion species (see the figure legend), solid and 
dotted lines represent spectra obtained on TPl and TP2 respectively, (c) Ion spectra obtained on TPl from 50 nm Cu target 
irradiated by a CP pulse at io ~ 1.25 x 10 20 W cm -2 . The proton spectra obtained with RCF, in the same shot, at different 
positions in the beam, corresponding to <~ 2°, ^7° and ~ 13° from the laser axis are also shown (dashed lines). 



quarter wave plate placed between the focussing optics 
and the PM. Here, e is the ratio between the laser electric 
field amplitudes along the vertical and horizontal axes, 
as determined by the combined effect of the wave plate 
and of reflection from the PM. Targets of different com- 
position and thickness were irradiated. Ion spectra pro- 
duced by the interaction were diagnosed by two thomp- 
son parabola (TP) spectrometers (with acceptance cone 
of 19.8 nSr) as shown in the Fig.[]Ja) - TPl placed along 
the laser axis and TP2 at 13 d= 2 degrees off axis. The 
energy-resolved spatial profile of the bottom half of the 
ion beam was recorded by stacks of radiochromic films 
(RCF). The image plate detectors used in the TPs were 
cross-calibrated with CR39 solid state nuclear track de- 
tector for particle number density [2l|. Similarly, the 



RCF dose response was absolutely calibrated [22]. 

While exponential spectra were always produced from 
5-10 fim thick foil targets, narrow-band features in pro- 
ton and heavier ion spectra were observed from sub-/im 
thick targets irradiated with high intensities. For exam- 
ple, the spectra in Fig.[]Jb), obtained from a 100 nm thick 
Cu target irradiated by a LP (p-polarization) laser pulse 
at peak intensity of ~ 3 x 10 20 W cm -2 , shows narrow- 
band peaked features in the proton (charge to mass ra- 
tio Z/A = 1) and carbon (Z/A = 0.5) spectra clearly 
separated from a lower energy component (as usual in 
standard interaction conditions, hydrogen and Carbon 
ions observed in the spectrum originate from surface con- 
taminant layers). The lower Z/A ion species had broad, 
exponential-like spectra similar to that shown in Fig.l 
(c) for Z/A ~ 0.42. Such charge state is likely to cor- 
respond to partially ionized Cu, as, for instance, a close 
inspection suggests a fine structure within the track pos- 
sibly corresponding to closely spaced values of Z/A. The 
peaks in ion spectra appear to be ordered, with the pro- 
ton peak at slightly higher energy than the C peak. It 
is also interesting to note the ion flux at the peaked fea- 



tures, which is nearly an order of magnitude higher than 
previously reported 161-18]. 

The narrow-band spectral features observed in carbon 
spectra along the laser axis are not observed with the 
off axis detector TP2, as shown in Fig. QJb), indicating 
that they are confined within a cone of half-aperture less 
than 13 ±2 deg. This indication is corroborated by RCF 
stack data taken simultaneously to the spectral measure- 
ment with the TPs (see Fig. (He)). The stopping range 
of carbon is significantly shorter than for protons at the 
same energy/nucleon: for example, the stopping range 
of 5 MeV/nucleon proton and carbon in the RCF stack, 
protected by a 30 /im Al foil, are ~ 240/im and ~ 60/im 
respectively. For this reason, the deposited dose in the 
RCFs in the stack can be assumed to be primarily due 
to the protons. By deconvolving the proton spectrum 
from the RCFs, the narrow band feature in the proton 
spectrum observed in the TPl was reproduced (dashed 
line in Fig. G^c)) for a defined region in the RCF corre- 
sponding to a half cone beam divergence of ~ 10°. The 
proton spectrum gradually becomes exponential (filling 
the dip in the spectrum at the low energy side) as one 
moves farther from the laser axis. Such exponential spec- 
tra are typical of TNSA acceleration caused by the ex- 
tended sheath produced by hot electrons at the target 
rear surface Assuming that the divergence of the 
narrow band feature in the proton beam and in the car- 
bon ion beam are comparable, the conversion efficiency 
into this component can be estimated as ~ 1%, which is 
significantly higher than reported in [18| and comparable 
to ref.0. 

The appearance and position of distinct peaks in the 
ion spectrum could be controlled by varying laser and 
target parameters as shown in Fig. [2(a) and (b). Peaks 
were only observed in the limit of thin foils and high 
intensity (Fig. 12(a)) with the peaks shifting towards high 
energy as either the intensity was increased or the target 
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FIG. 2: (a) Graph showing comparison between five ion (Z/A 
= 0.5) spectra, where the spectral peaks are plotted in (b) 
as a function of a^r p /x- The experimental parameter set [ao, 
target material, target thickness (/im)] for the data points 1-7 
are [15.5, Cu, 0.1], [10, Cu, 0.05], [13.8, Cu, 0.1], [7.5, Al, 0.1], 
[6.9, Al, 0.1], [13.6, Al, 0.5] and [14.1, Al, 0.8] respectively. 
Black solid line in (b) represents ion energy estimated from 
the analytical modeling described in the text. The red circle 
represent the data shown by Henig et. al. in Ref. [l6[. 



thickness was reduced (Fig. EJb)) 

In order to assess the possible influence of Radiation 
Pressure effects on the spectral profiles observed, a sim- 
ple analytical model was developed taking into account 
the Hole Boring (HB) and LS phases of RPA mecha- 
nism. Due to the extreme radiation pressure (of the or- 
der of hundreds of GBar) exerted on the irradiated target 
surface, ions are swept forward by directional momen- 
tum transfer d-Q- Consequently, the laser pulse bores 
through the target, in the HB phase, with the ion front 
velocity (in units of c) given by [H, 0, [HI 
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where a(t) = 0.85a// (t)A 2 /10 18 W cm" 2 , p' = p/m p n c 
with p the mass density, m p the proton mass and n c is 
the critical electron density, and Io(t) is the laser inten- 
sity on as a function of time. (Relativistic effects Q are 
neglected here for simplicity, consistently with the values 
of j3 <C 1 for our parameters.) If the target is sufficiently 
thin, the hole boring front will tend to reach the target 
rear surface at a time tuB < t p (laser pulse duration), 
i.e. before the end of the laser pulse. As the thickness 
of the compressed layer becomes comparable or less than 
the evanescence length of the ponderomotive force, the 
whole layer can be cyclically accelerated with high effi- 
ciency for the remainder of the laser pulse duration 
In this scenario of whole foil acceleration (LS regime) the 
equations of motion of the foil can be expressed as [8l4l2j 
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where r = ct/X, 7 = (1 — f3 2 ) l / 2 and x = p'l/X for a tar- 
get of thickness /, under the approximation that the areal 



density of the compressed layer is same as that of the tar- 
get before the interaction. The reflectivity R = R((3, x> °) 
may be estimated using a simple model of relativistically 
induced self-transparency ^ and the system of Eqs.(0) 
can be integrated numerically. However, some insight 
may already be obtained analytically. Assuming R ~ 1, 
Eq.(j2]) may be integrated to yield the final velocity of 
the foil /3f and the corresponding energy per nucleon 
Eion = (7/ — l)m p c 2 as a function of the dimension- 
less fluence parameter £ = (pic 2 ) -1 Ip(t)dt and of 
the initial velocity f3(t ) = /3 Q: 
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To match the initial HB stage with the later LS one, we 
take to = tnB and, from Eq.(pQ), /?o = yHtoJ/pc 3 - For 
thin targets and high intensities, tjjB < ^ —^ls and the 
LS stage will dominate also because of higher efficiency 
with respect to HB, thus £ ^> Zq. In such limit we may 
estimate £ ~ OqTls/x — a o r p/x where r p = ct p /X and 
a o r P ~ I^° 00 a ' 2 ( r ) dr, and E ion will scale as (do r p/x) a 
with a = 2 for £ <C 1 and a = 1 for £ ^> 1, i.e. in the 
ultra-relativistic case jjjj. 

On the basis of the above equations, we can estimate 
whether the peaks observed, and their scaling with the 
target and laser parameters are consistent with the ex- 
pectations for LS. The switch-on time for the LS phase 
is taken as the time at which the compressed ion front 
reaches the target rear surface as a result of the HB pro- 
cess. In this simple rigid model, the degree of ionisation 
of target ions and the number of different ion species 
present in the compressed layer are irrelevant. There- 
fore, without any loss of generality, we have ignored the 
target contaminant layers, which have typically few nm 
thickness and significantly lower density than the metal 
targets used in our case. As shown in the Fig. |2jb), the 
ion energy estimated by the model agrees well with the 
data, within the experimental errors. As expected for 
the non-relativistic case, the ion energy is seen to scale 
as (doTp/x) 2 - For comparison, the only published data 
point for carbon peaks attributed to RPA-LS [l6[ is also 
shown in the graph, which is in substantial agreement 
with the calculated scaling. For the data points shown 
in the Fig. [21(b), the hole boring phase ends significantly 
before the peak of the laser pulse and the non-linear re- 
flectivity remains at unity for the duration of the pulse. 

The experimental features described above are ob- 
served virtually independent of the incident laser polar- 
ization. This is not surprising since the cycle-averaged 
light pressure is the same for CP and LP. While the early 
simulations of RPA were carried out employing LP pulses 
at extreme intensities [H, E|; the attention of theoreti- 
cal work over the past few years has focused on the use 
of circular polarization as a means to minimize effects 
of electron heating and sheath acceleration, and achieve 



a stable LS drive. Multi-dimensional simulations show 
that the stability of the LS acceleration is dominated by 
factors such as the shape of the radial laser intensity pro- 
file [sl, [l3j and achieving a smooth transition from the HB 
phase to the LS phase at the early stages of the acceler- 
ation process [25]. The polarization, while still affecting 
the outcome, plays a sub-ordinate role in this case. One 
should note that, in the experiment the effect of chang- 
ing the polarization of the laser pulse may be weakened 
by several factors, such as the polarisation-dependent re- 
flectivity of the plasma mirror (so that circular polariza- 
tion actually corresponds to an elliptical polarisation on 
target with e ~ 1.14), the tight f/3 focusing resulting 
a spatially Gaussian intensity profile on target, and by 
the relatively long duration of the laser pulse. The lat- 
ter two factors lead to a significant deformation of the 
ultra-thin target foil during the acceleration so that nor- 
mal incidence is not strictly preserved locally. All these 
factors enhance electron heating and reduce any polar- 
isation dependance. Furthermore, recent work by Qiao 
et al. [25] has highlighted that efficient LS drive can be 
achieved with LP pulses at presently achievable inten- 
sities, in a regime where RPA and TNSA co-exist. In 
this framework, sheath-field effects may account for the 
observed separation of the C and proton spectral peaks 
[see Fig. 0Jb),(c)] (as in an ideal LS scenario all ions in 
the sail should reach the same energy per nucleon). Af- 
ter the LS stage, the sheath field can in principle further 
accelerate the ions, with the protons gaining more en- 
ergy /nucleon than Carbon ions due to their higher Z/A 
. This last acceleration stage may preserve or even en- 
force the peaked distribution as it typically happens in a 
multispecies expansion [26|, and references therein]. 

This scenario is confirmed by 2D PIC simulation em- 
ploying the "ILLUMINATION" code ^ . Due to limita- 
tion in computational resources which prevent ~ps scale 
simulation of long pulse interaction, we scaled down the 
laser and target parameters in such a way that both the 
dimensionless quantities Io/pc 3 and («o r p/x) remain the 
same as the case shown in the Fig.QJb). A LP laser pulse 
with A = 1.0/im, I = 7.5 x 10 19 W cm" 2 (a = 7.75), 
gaussian spatial profile of 5/im radius and gaussian tem- 
poral profile of 50 fs FWHM was used. A copper target 
was used with electron density of 384. 3n c and thickness 
of Zo = 10 nm, covered on both side with 3 nm, 48.0n c 
layer of hydrocarbon contaminants. The simulation box 
was 20 /im long and 24.576 /xm wide with 1 nm x 8 nm 
cell size in order to assure resolution for the contaminant 
layers. 72 particles per cell for Cu 28+ species and 36 par- 
ticles per cell for H + and C 6+ in the contaminant layers 
were used. 

As expected, the simulation shows that TNSA and HB 
takes place with the rising intensity of the laser pulse, ac- 
celerating ions from rear and front surfaces respectively. 
After the ions from the front surface pile up at the tar- 
get rear surface forming a compressed layer, ions are ac- 
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FIG. 3: (a) Ion spectra obtained from 2D multilayered PIC 
simulation for the data point '1' in Fig. [2ja), of which the 
experimental spectra is shown in the Fig. [TJb). (b) and (c) 
shows 2D density profiles of ions (black-proton, red-carbon 
and blue- Copper) at 66 fs and 165 fs, respectively. Each col- 
ormap ranges from zero to the value mentioned in the graph. 



celerated predominantly by the LS mechanism resulting 
in the ion spectrum shown in the Fig. [3ja), which re- 
produces the main features observed in the experimental 
data. The target remains highly reflective during the 
whole duration of the simulation and the carbon ions are 
accelerated in a snowplough fashion maintaining a dense 
layer as shown in the Fig. [3]^b) and (c). It is interest- 
ing to note the debunching of the Cu ions after the laser 
pulse [12], as shown in the Fig. [3fc), producing the ex- 
ponential spectrum profile of Fig. [3fa), which resembles 
closely the experimental data. This is consistent with the 
scenario highlighted in [11], where the heavier species in a 
multispecies target undergo decompression and Coulomb 
explosion, while the lighter species are stabilized by the 
excess electrons released by the heavier species. 

The experimental data points in Fig. Efb) show that 
the ion energy scales with the parameter («o t p/x) to the 
power 2, as expected for a fully reflecting target, and 
in the limit /3 <C 1. As reported by several groups via 
extensive 2D and 3D simulations (see Fig. 2]), the ion 
energy can therefore be enhanced by increasing laser flu- 
ence and decreasing the target areal density. However, in 
doing this, one needs to avoid self-induced transparency, 
as this terminates efficient LS (see the blue lines in the 
Fig. 2]) and leads to a dramatical reduction of the ion 
energy. 

Although the ion energies achieved here (>10 
MeV/nucleon) are encouraging, producing peaks at more 
than 100 MeV/nucleon would be a crucial milestone in 
view of applications. As shown in the insert in the Fig.HJ 
2D PIC simulations predict that 100 MeV/nucleon ions 
in a narrow energy bandwidth can be reached by, for ex- 
ample, increasing the laser fiuence by a factor of 2 and 
decreasing the target density by a factor of 2.5, compared 
to the case shown in the Fig.[TJb). This seems achievable 
given current developments in laser technology. 

In conclusion, we have reported on the observation of 
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FIG. 4: Ion energy scaling (black solid line) as shown in the 
Fig. [2jb), extrapolated to higher a§r p /x, assuming R = 1. 
Square points are the experimental data shown in the same 
figure. The red circles represent the results reported in the 
literature, as labeled, by multispecies PIC simulations for sta- 
ble LS acceleration. Dashed black line shows trend for a 
Eion oc (a>oTp/x) 2 scaling valid for non-relativistic ion ener- 
gies. Solid and dashed blue lines plots expected ion energy 
(according to the rigid model) by varying x f° r 45 fs FWHM 
laser at J = 5 x 10 20 W cm" 2 and 450 fs FWHM laser at 
io = 5 x 10 19 W cm -2 respectively. The insert shows the 
spectra obtained from a PIC simulation identical to the run 
shown in the Fig. [2jc), except the laser fluence increased by 
a factor of 2 and target density decreased by a factor of 2.5. 

narrow band features in the spectra of laser-accelerated 
ions, which appear to be consistent with Radiation pres- 
sure acceleration, in a regime where LS overcomes sheath 
acceleration. The observed ion beams show unique prop- 
erties in terms of ion energy, divergence and fast energy 
scaling, and offer high promise for further progress in the 
near future. 
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